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Abstract. The economic life of an asset is the optimum length of its use- 
fulness, which is the moment that the asset's expenses are minimum. In this 
paper, the economic life of physical assets, such as industry machine and equip- 
ment, can be interpreted as the moment that the minimum is reached by its 
equivalent property cost function, defined as the sum of all equivalent capital 
and maintenance costs during its life. 

Many authors in classical papers have used principles of engineering eco- 
nomic to solve the assets replacement problem. However, in the literature, the 
main attributes found were proved with intuitive ideas instead mathematical 
analysis. Therefore, in this paper the main goal is to study these principles of 
engineering economic with mathematical techniques. 

Here, is used non-smooth analysis to classify all the possibilities for the 
minimum of a class of equivalent property cost functions of assets. The mini- 
mum of these function gives the optimum moment for the asset to be replaced, 
i.e., its economic life. 



1. Introduction to physical assets replacement 

Physical assets, such as industry equipment, are vulnerable to devaluation and 
obsolescence. Among the many consequences of devaluation, there are, for instance, 
output decline, operation and maintenance expenses increase. These causes can lead 
to market value declination (see Park and Sharpe-Bette in T)- Regarding asset's 
obsolescence, the main causes are technological innovation and change inside the 
company organization. Therefore, the asset is no longer needed (see Park in [B]). 

Nevertheless, with the proper maintenance, physical assets can be used for much 
more time than its physical nature allow. For example, it is possible to see vintage 
cars driving on the streets. However, to make it possible, the companies must be 
willing to pay a higher price. 

Along with that comes the concept of Economic Life. An asset's economic life 
is the length of its usefulness, in a way that the expenses, i.e., annually sum of the 
maintenance costs and capital costs, are minimum. Therefore, the economic life of 
the asset is the optimum moment to replace the asset. 

Because of that, if the asset is kept longer than its economic life, the expenses of 
maintenance will have increase a lot. Meanwhile, if the asset is replaced before its 
economic life, the capital cost will not have been fully fiscal depreciated. Therefore, 
part of the investment, in the acquisition cost of the asset will be lost. So, physical 
assets in general are always used for a limited time. 



Key words and phrases, physical assets, replacement problem, economic life, non-smooth 
analysis. 
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2. Objectives 

Many authors, in classical papers, like Alchian in ^J; Park in Park and Sharp- 
Bette in [7]; Grant, et al in and Thuesen, et al in [5], have used principles of 
engineering economic, as net present value and annuity equivalent, in the assets 
replacement problem. 

However, in the literature, the main attributes found were proved with intuitive 
ideas and with no mathematical accuracy. Therefore, in this paper the main goal 
is to study these principles of engineering economic with mathematical analysis. 

3. Mathematical modeling 
The concept of economic life can be modeled with the following equations: 

3.1. Equivalent Capital CosiQ. Two components arc involved in the evaluation 
of capital cost: Acquisition Cost A > 0, and Salvage Value R(t). The first is a 
fixed value, while the second changes along the time. The result obtained from 
the difference between the annual equivalent and each component is the equivalent 
capital cost (see Appendix A), as showed in equation ((l}. 

(1) ECC g{t) = 1^ {Ae^' - R{t)) . 

Here, < r < 1 is the nominal interest rate. It is usual to see the value around 
10%. 

3.2. Equivalent Maintenance Cost0. To compute the equivalent maintenance 
cost, firstly it is necessary to evaluate the sum of present value of the series of Main- 
tenance Costs M{t) along the time. Then compute the present value in annually 
equivalents, as it is showed in equation ([2]). 

(2) EMC fit) = ^^e^' f M{s)e-^^ ds. 

e — ^ Jo 

3.3. Equivalent Property Cost. The equivalent property cost of an asset is the 
sum of all capital and maintenance equivalent costs during its economic life, as it 
can be seen in equation ([3|). 

(3) EPC''^h(t)=f{t)+g{t). 

It is possible to see in Figure [1] the behavior of the functions ([1]), ©, and ([3]) 
along the time. 

The Figure [1] is representing the costs (y-axis) of the functions, ([H, ([2]) and 
(|ni), along the time (x-axis). It can be seen that the equivalent maintenance costs 
increase, while the equivalent capital cost decrease. Therefore, the sum of both 
variables, i.e., the equivalent property cost, reaches a minimum point, which is the 

^The Capital Cost Acquisition is the cash outflow regarding asset's purchase. This cost also 
include shipping cost, installation cost and training cost. Besides, this is also called first cost or 
investment cost, because this kind of expense is the cost of getting an activity or project started 
(see Fabrycky and Blanchard in [3], p. 22). 

■^On this category the main costs are due to cleaning, lubrication, components adjustment and 
repair, among others (see Hastings in [5]). 
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Figure 1. Model to estimate the economic life. 

asset's economic life. In Figure [1] the minimum point is reached at year five. So, 
it is the optimum moment for the asset be kept. 

On balance, the consequences of a late replacement imply in high expenditures 
of operational and maintenance costs. Besides, there is also the opportunity cost 
of loss of asset's market value. On the other hand, an early replacement imply in 
selling the asset before its capital recovery . 

4. Physical Asset Model 

The maintenance costs function M{t) is usually taken as a linear regression from 
a data collect regarding maintenance costs within the company. Therefore, in many 
situations, it can be accepted that M{t) is linear increasing function (see (jH)), since 
the maintenance cost are higher each year. 



* ^ -A-EPC 




i-EMC 



•♦■ECC 



(4) M{t) = at, 

with a > 0. 

The salvage value function R{t) is the market value of the asset. It can be 
expected that R(t) is decreasing, since physical assets are vulnerable to devaluation. 
However, the market value is always positive. So, given A, the acquisition cost 
of the asset, it is possible to assume, in many situations, that the asset value is 
depreciated every year by a tax 6, i.e., R{t) — A — bt, until its fully depreciate at 
t = A/b. Thereafter, the asset value is constant equal zero (see ([5])). 



(5) R{t) = 

with A> and 6 > 0. 



A-bt if < i < A/b, 
if t > A/b, 
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With the assumptions above, the equivalent property cost function ([3]) becomes 
(e*^ - 1) frt{br - a) 



(6) hit) 



6"^* - 1 



+ a + Ar-'] if < t < A/b, 



^l^l±^(e-(a + Ar^)-a(l + rt)) if t > A/b^ 



and its derivative, for t ^ A/b, is given by 
(e'- - 1) 



(7) h'{t) = 



) r(e--lf 



{e''\rt - 1) + l) {a - br) if 0<t<A/b, 



{a - e''' {-art + a + Ar'^)) if t > A/b. 



r(e'-* - 1) 

5. Statements of the main results 
For u > 0, define the function 

r(M) = 1 + li + Wo (-e-i-") , 

where Wq represents the main branch of the multivalued Lambert W-Function (see 
Appendix B). 

Our main results, that classify all the possibilities for the minimum of the func- 
tion h{t), are the following: 

Theorem 1. If c = Ar'^ / a, then we have the following possibilities for the minimum 
of the function h{t): 

(Ci) z/ 

Abr'^ 



a> br and a > 



Ar + b (^e^'ir - ij 

then the minimum of the function h(t), for t > 0, is uniquely reached at 
t = 0. 

{C2) ^/ 

Abr"^ 



a = br and a > 



Ar + b{e 5^-1 



then the minimum of the function h{t), for t > 0, is reached at all t G 
[0,A/b]. 

(C3) tf 

Abr^ 

a < br and a> — 



Ar + b\ 



then the minimum of the function h(t), for t > 0, is uniquely reached at 
t = A/b. 

(C4) if 

Abr^ 



a> br and a < 



Ar + b (^e'^ - ij 



then the minimum of the function h(t), for t > 0, is reached at t — 0, or 
t = T{c)/r. 
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a < br and a < 



Abr^ 



Ar + b[e r - l 



then the minimum of the function h(t), for t > 0, is uniquely reached at 
t = T{c)/r. 



The proof of Theorem [T] is given in Section 7. The next theorem divides the case 

4 , C4 and C4 . 



C4, of Theorem [U in three different cases, C|, C| and 



Theorem 2. Considering the assumptions of the case C4 of Theorem\^ follows: 

{CD If 

A > log ( 1 - - ) - 

then the minimum of the function h(t), for t > 0, is uniquely reached at 
t = 0; 



ici) If 



A = log (1 - - 
a 



then the minimum of the function h(t), for t > 0, is reached at t — and 
t = T{c)/r; 



{CD ^f 



A < log ( 1 - -, 

a I r 



then the minimum of the function h(t), for t > 0, is uniquely reached at 
t = T{c)/r. 

The proof of Theorem [2] is given in Section 7. 

Resuming the results of Theorems [T] and [21 it follows the classification: 

• in the cases Ci and C\ (see Figure [5]), the minimum is uniquely reached at 
t = 0; 





A/b T{c)/r 



Figure 2 . Cases Ci and C] . 




Figure 3. Case C2. 

in the case C2 (see Figure [3l), the minimum is reached at ah t e [0 , 

in the case C3 (see Figure|4]), the minimum is uniquely reached a.t t = A/b; 
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Figure 4. Case C3. 

in the case C4 (see Figure [5]), the minimum is reached at two points t — 

and t = T{c)/r; 

h{t) 
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Figure 5. 


Case C|. 
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in the cases Cf and C5 (see Figure [6]), the minimum is uniquely reached at 

t = r(c)/r; 
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Figure 6. Cases C4 and C5. 



6. Basic lemmas 

Before proving the Theorems [1] and [H we shall prove two lemmas about the 
critical points of the function © on the intervals Ii = (0, A/b) and I2 — {A/b, +00). 
For this aim, we study the zeros of the derivative ([7]) . Observe that the functions 
h{t) and h'{t), given respectively in ([6]) and ([7]), restricted to the intervals Ii and 
I2, are smooth. 

Lemma 3. If Ii — {0,A/b), then h'{t) ^ for t E Ii. Moreover, if a > br, then 
h{t) is increasing in Ii; if a < br, then h(t) is decreasing in Ii; and if a — br, then 
h{t) is constant in Ii. 



Proof. For t e Ii 



(e'' - I) 



r(e''* - 1)" 

where Pi{t) = (a - &r) (1 + e^(r - 1)). Thus for t E Ii, we have that h'{t) = if 
and only if pi(rt) = 0. 

On the other hand, for a ^ br, Pi{t) = if and only if t = 0. Therefore, for 
t E Ii and br ^ a, h'{t) = if and only if < = 0. 

Now, observe that, for t E h, 

> 0, if a > br, 
h'{t) { =0, if a = br, 
< 0, if a < br. 



Moreover, if a = br, then, for t E h. 



hit) 



e'' - 1 



(6 + Ar) 



ioT < t < A/b, which concludes the proof. 



□ 



Lemma 4. For t > A/b follows: 
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(i) ^I 

Abr^ 

Oj > J r-, 

Ar + b(e-^ - 1 j 
then h'{t) ^ 0, for t > A/b; Moreover, h{t) is increasing for t > A/b; 

(ii) ^f 

Abr"^ 

A + b{e-T - ij 

then h'{t) =0 fort = T{c)/r, where 

l + c + Wo{-e-'-'), 

where c = Ar'^/a. Moreover, h{t) is decreasing for A/b < t < T{c); and 
h{t) is increasing for t > T{c); 

Proof For t > A/b, 

r (e'^* — 1) 

where P2{t) = a — e'^{a + Ar'^ — ar). Thus for t > A/b, we have that h'{t) = if 
and only if P2{rt) =0. 

On the other hand, P2{rt) = if and only if p{rt) = c, where 

re'^ -e'^ + 1 , Ar"^ 
P{t) = and c = . 

e"^ a 

Note that the function p(t) is increasing for r > 0. Indeed p'{t) = 1 — cosh(T) + 
sinh(T). Thus if p(rto) = c, then tg > A/b if and only if 

/ A\ Ar"^ Ar Abr'^ 
p{rto)>p\r—\ >— +e "S" - 1 <s=> a< 



A + b{e-T 



Hence if 



then h'{t) ^ Q ior t > A/b. 
Now, assume that 



Abr"^ 

a > 



Ar + b{e~^ - 1 j 



Abr"^ 

a < 



Ar + b{e-^ -l)" 
Using the software Mathematica 8.0, we obtain that: if 



T(w) = 1 + u + Wo (- 



-e 



-l-u\ 



where Wq represents the main branch of the multivalued Lambert W-Function (see 
Appendix B), then 

T(c) 



c. 

r 

Again, using the software Mathematica 8.0, we can check that T{u) > for 
u > 0, which implies, as we have seen, that T{c)/r > A/b. 

Now, using the fact that the function p(r) is increasing for t > 0, we have, in this 
case, that t = T{c)/r is the unique critical point of the function h(t) for t > A/b. 
Moreover, is easy to see that p{r t) < cfoi t < T{c) / r, and p{r t) > ciov t > T{c) / r. 
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Thus the derivative h'{t) is negative for A/b < t < T{c)/r, and is positive for 
t > T{c)/r. Hence we can conclude that h{t) is decreasing for A/b < t < T{c)/r; 
and h{t) is increasing for t > T{c)/r; which concludes the proof. □ 

7. Proofs of Theorems [U and [2] 

Now, we are ready to prove the Theorems [1] and [2] 

Proof of Theorem QJ The function (j6|) is continuous for t > and piecewise differ- 
entiable. Thus we can find the minimum of ([6|) comparing the minimums of each 
piecewise. In other words, we have to choose the minimum among the values of 
the function evaluated at the critical points of each piece and at the extreme points 
< = and t = A/b. 

In the Case Ci, by Lemmas [3] and IH the function h{t) is increasing for < > 0. 
Hence the minimum is uniquely reached at t = 0. 

In the Case C2, by Lemma [3l the function h{t) is constant for < t << A/b; 
and by Lemma [H the function h{t) is increasing for t > A/b. Hence the minimum 
is reached at all t G [0,v4/6]. 

In the Case C3, by Lemma[3l the function h{t) is decreasing for < t << A/b; 
and by Lemma [H the function h(t) is increasing for t < A/b. Hence the minimum 
is uniquely reached at t = A/b. 

In the Case C4, by Lemma [31 the function h{t) is increasing for < i < A/b, 
thus the minimum for < t < A/b is reached at i = 0. On the other hand, by 
Lemma S] there is a local minimum in t = T{c)/r > A/b. Moreover, by Lemma SI 
the function h(t) is increasing for t > T(c)/r. Hence the minimum is reached at 
t = 0, or t = T{c)/r. 

In the Case C5, by Lemma [3l the function h{t) is decreasing for < t < A/b, 
and by Lemma [31 there is a local minimum in t ~ T[c)/r > A/b. Moreover, by 
Lemma [31 the function h[t) is increasing for t > T{c)/r. Hence the minimum is 
uniquely reached aX t = T{c)/r. 

□ 

Proof of Theorem [H Consider the assumptions of the case C\ . It is clear that 
the minimum of the function h{t), for t > 0, is reached at t = if and only if 
/i(0) < h{T{c)/r). We know that 

h{Q)=^-—^{Ar + b) and /i (^^^ = (a + Ar^ + aWo (-e"^-"^)) . 

Thus /i(0) < h{T{c)/r) if and only if 

(8) _ 1 < ^ _ 1 < T4.„ (_e-i-^) . 

The function Wo : (—1/e, 00) — > (— l,oo) is invertible and Wq^ : (— l,oo) — !■ 
(—1/e, 00) is an increasing function (see Appendix B). Then we can apply the 
function Wq^ in both sides of the inequality ([5|), since br/a — \ > —1, to obtain 
the equivalent inequality 



(9) 



1 - 



br 

br\ — -1 

e a 



> e 



> 0. 
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Now, applying the log function in both sides of the inequality ([9]), we can con- 
clude that h{0) < h[T{c)/r) if and only if 

. , A 6r\"^ h 

Hence we conclude the proof of theorem for the case C\ . The proofs for the cases 
C| and C|, are completely analogous. 

□ 

8. Conclusion and future directions 

The methodology defined in the classical literature, assets replacement prob- 
lem involve solving it organizing and analyzing tables and equations, as it can be 
verified. However, this involve too much work and effort. So, for assets with the 
characteristics described in Section 2] it can be applied the methodology described 
in this paper. With the analytical approach, it takes less time and effort to see 
the when is the economic life and to replace the asset. Following this article ap- 
proach, through Theorem [1] and Theorem [2] results, it can simply solve the assets 
replacement problem. 

The function M{t) in the expression ([2]) comes from a linear regression. Thus 
in a more general case, M{t) could be a piecewise function composed by affine and 
exponential translated functions. For instance (see Figure [7]), we can take 

(10) Af(t) = ^X7.(i)M^(t), 

where, iV is a collection of natural indexes; is a sequence of disjoint real 

intervals such that covers the real line M; for each i N, Mi is an affine, or 
exponential translated function; and x/ : R — ?> {0, 1} is the characteristic function, 
i.e., if / is a subset of M, then 

Xiit) = 

We intend, in future works, deal with the problem of assets replacement consid- 
ering maintenance costs functions as given in ^TU\\ . 

Appendix A: Financial Criterion Equivalence 

Given a series of irregular cash flow from i = 1 to A , we want to make it a 
regular cash flow along the years. The scenario is possible to see in Figure [8] below. 

In Figure IH it is possible to see that first we have a cash flow with different 
values along the years. Then, we have a new cash flow with equal deposits along 
the years. Now, the question is how to make this equivalence? 

Let us consider a future deposit F at time A and a series of regular deposit A 
from year 1 to A. Then, to make the equivalence, it is necessary to consider the 
interest i for each year. Therefore, 

N 

(11) F = A^(l-f ^)^-^ 
Multiplying equation (fTTj) by (1 + i), follows 
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M{t) 




Figure 7. Piecewise function M{t). 
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Figure 8. Equivalence of cash flows. 



(12) 



N 



+ N-k 



k=l 



Now, computing the difference between equations p2)) and (|lll) . we have that 

iF = A{{l + if-l). 
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Since F — P{1 + ?)^, we can conclude that 



(13) '—jrf A- , 

As a consequence of equivalence (|13p , we can take any present value and trans- 
forming in a series of regular cash flows. 

Consider r the nominal interest, M the number of period per year, and i the 
real interest. It is possible to establish the equivalence 



(14) .= (l + ^)"-l. 

Computing the limit of the expression (|14p . when M goes to +oo, it follows that 
i = eJ — 1. Hence, it can be establish the equivalences between two discount factors: 

e-- = ^— 
(1 + z)*- 



Appendix B: Lambert W Function 

Let ^ G C be any complex number and define VF(z) C C as a set of complex 
numbers w, such that satisfies the equation 

w 

z = we . 

W(z) is called Lambert W Function. It is defined to be the multivalued inverse 
of the function w M> we™, for w £ C 

If z is real, then for— l/e>z<0 there are two real numbers in the set W{z) 
(see Figureini). We denote the branch satisfying —1 < W{z) by Wo{z), or just W{z) 
when there is no possibility for confusion, and the branch satisfying W{z) > — 1 by 
W-iiz). 

Wo (2) is referred to as the principal branch of the W function, and it is analytic 
at 0. Moreover, the series expansion for Wo{z) is given by 

(15) Wo{z)^Y.-^'"^ 

n=l 

which has the radius of convergence equal to 1/e, i.e., the series converges for 
ah z e (-1/e, 1/e). 

If we consider the restricted function Wq : (—1/e, 00) — > (— l,oo), then Wq 
becomes invertible with Wq~^ : (— l,oo) — > (—1/e, 00) being an increasing function. 
Indeed, Wq^{w) — we™, for w G (— l,oo). 

For more detail on the Lambert W function, see, for instance, [2]. 
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Figure 9. The two real branches for W{z). , Wo{z); , W-i{z). 
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